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We study an optimum measurement for quantum state discrimination, which maximizes the prob- 
ability of correct results when the probability of inconclusive results is fixed at a given value. The 
measurement describes minimum-error discrimination if this value is zero, while under certain condi- 
tions it corresponds to optimized maximum-confidence discrimination, or to optimum unambiguous 
discrimination, respectively, when the fixed value reaches a definite minimum. Using operator condi- 
tions that determine the optimum measurement, we derive analytical solutions for the discrimination 
of two mixed qubit states, including the case of two pure states occurring with arbitrary prior prob- 
abilities, and for the discrimination of N symmetric states, both pure and mixed. We also consider 
a case where the given density operators resolve the identity operator, and we specify the optimal- 
ity conditions for partially symmetric states. Moreover, we show that from the complete solution 
for arbitrary values of the fixed rate of inconclusive results one can always obtain the optimum 
measurement in another strategy where the error rate is fixed, and vice versa. 
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I. INTRODUCTION 

In quantum state discrimination we want to deter- 
mine the actual state of a quantum system that is 
known to be in a certain state belonging to a given 
set of states. This is an essential problem for many 
tasks in quantum communication and quantum cryptog- 
raphy. Since nonorthogonal quantum states cannot be 
distinguished perfectly, various optimized discrimination 
strategies have been developed. The best known of these 
are minimum-error discrimination 0, Q and optimum un- 
ambiguous discrimination, where the latter strategy was 
originally derived for pure states 0-0] and later also con- 
sidered for mixed states In unambiguous discrim- 
ination errors do not occur, that is, the total error rate 
P e is required to vanish, P e — 0. This can be achieved at 
the expense of admitting inconclusive results, where the 
measurement fails to give a definite answer. Optimum 
unambiguous discrimination minimizes this failure prob- 
ability, or failure rate, Q, yielding the minimum value 

Q Qmin • 

Unambiguous state discrimination is not possible for 
pure states that are linearly dependent [l5| , or for mixed 
states described by density operators with identical sup- 
ports 0, H|, where the support is the Hilbert space 
spanned by the eigenvectors with non-zero eigenvalues. 
When unambiguous discrimination is impossible, related 
measurements can be applied, which discriminate the 
states with maximum possible confidence for each conclu- 
sive outcome (la . Optimized maximum-confidence 
discrimination is achieved by the particular one of these 
measurements that minimizes the failure probability [l8l - 
|22T |. This measurement corresponds to optimum unam- 
biguous discrimination when for each of the different con- 
clusive results the maximum confidence is equal to unity. 

While the confidence is defined separately for each con- 
clusive outcome, in minimum-error discrimination the 



overall error rate P e , averaged over all outcomes, is mini- 
mized. In this strategy inconclusive results are not al- 
lowed, that is, Q = 0. Chcfles and Barnett [23[ in- 
troduced a more general strategy that minimizes P e 
when a certain fixed rate Q of inconclusive results is ad- 
mitted, thus reducing the minimum achievable value of 
P e . For discriminating two equiprobable nonorthogonal 
pure states these authors obtained the optimum solution, 
which interpolates between minimum-error discrimina- 
tion and optimum unambiguous discrimination when the 
fixed rate Q grows from Q = to Q = Q m in- Their re- 
search triggered further investigations |24| - l29| . In partic- 
ular, optimum state discrimination with a fixed rate Q of 
inconclusive results was extended to mixed states [2fj| . [26| . 
and the relation to the strategy of maximum-confidence 
discrimination was briefly discussed in our previous pa- 
pers [HI, [22j • Clearly, in a measurement where Q is fixed, 
minimizing P e corresponds to maximizing the overall rate 
of correct results, P c . Fiurasek and Jezek 25| derived 
general operator conditions that have to be fulfilled in 
the optimum measurement with an arbitrary fixed value 
Q, holding for the discrimination of an arbitrary number 
TV of mixed states. For N — 2 they solved the optimiza- 
tion problem in the special case of two mixed qubit that 
have the same purity and occur with equal prior proba- 
bilities [25J , which includes the solution for two equiprob- 
able pure states [23|. Other solutions were not obtained. 

Recently optimum state discrimination has been also 
investigated for a measurement where a certain fixed er- 
ror rate P e is admitted [27i - [29| . With a fixed rate of P e , 
maximizing the rate of correct results P c corresponds to 
minimizing the failure rate Q. The optimum measure- 
ment in this strategy was studied for the discrimination 
of two states, N = 2, where the solution was derived for 
two pure states occurring with arbitrary prior probabil- 
ities [28|,|29|. However, it was not recognized that there 
exists a relation between optimum state discrimination 
with a fixed value of Q, on the one hand, and with a fixed 
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value of P e , on the other hand. This relation implies that 
by completely solving the optimization problem for one 
of these two discrimination strategies one can also obtain 
the solution for the other strategy, as will be shown in 
the present paper. 

The purpose of our paper is twofold. First, for opti- 
mum state discrimination with a fixed rate of Q we de- 
rive analytical solutions going beyond the solutions ob- 
tained so far [H, [25| , and we also outline the relation to 
maximum-confidence discrimination, extending our pre- 
vious discussions [H, Second, we show that when 
the maximum rate of correct results P c is known as a 
function of the fixed failure probability Q, then from this 
function one can also obtain the maximum rate of cor- 
rect results in dependence of a fixed error rate P e , and 
vice versa, due to a general relation between the solu- 
tions of the two optimized strategies that holds for the 
discrimination of an arbitrary number N of mixed states. 

The paper is organized as follows: We begin in Sec. 
II with an alternative derivation of the optimality con- 
ditions for state discrimination with fixed Q, and with 
establishing the relation to maximum-confidence discrim- 
ination. In Sec. Ill together with Appendix A we apply 
these conditions to the discrimination of two mixed qubit 
states occurring with arbitrary prior probabilities, while 
Sec. IV together with Appendix B is devoted to the dis- 
crimination of N symmetric states. In Sec. V we consider 
a case where the density operators resolve the identity op- 
erator, and we also specify the optimality conditions for 
a case of partially symmetric states. The general relation 
between the two optimization strategies where either P e 
or Q has a fixed value is presented in Sec. VI, together 
with an example. Sec. VII provides a summary of results 
and concludes the paper. 

II. OPTIMUM MEASUREMENT 
A. Conditions for optimality 

We suppose that a quantum system is prepared with 
the prior probability rjj in one of TV given states de- 
scribed by the density operators pj (j = 1, . . . , N), where 
Y]^—i T]j = 1. The general task is to perform a discrim- 
ination measurement in order to infer in which of the 
N possible states the system was prepared. The mea- 
surement is described by N + 1 positive detection oper- 
ators III, • • ■ IIjv and IIo, where Tr (ptJIj) is the condi- 
tional probability that a system is inferred to be in the 
state pj given it had been prepared in the state pk , while 
Tr (pfcllo) is the conditional probability that in this case 
an inconclusive result is obtained and the measurement 
fails to discriminate the states. The detection operators 
fulfill the completeness relation Hj = /, where I is 

the identity operator in the c£-dimensional Hilbert space 
Hd spanned by the eigenstates of the operators p%,..., pn 
that belong to nonzero eigenvalues. Unless each of the 
detection operators is a projector, the measurement is a 



generalized one. Once the detection operators are known, 
implementations of the generalized measurement as a 
projective measurement in an enlarged Hilbert space can 
be obtained using standard methods [30|, HH . 

We are interested in the specific discrimination mea- 
surement where the overall failure probability 

N N 

Q = ^2 ? 'j Tr (ft n o) = Tr G° n °) with P = ^2 ViPi 
i=i j'=i 

(2.1) 

has a given fixed value Q < 1, while the overall probabil- 
ity of getting correct results, P c , is as large as possible. 
Here we defined the total density operator p, which has 
its support in the full Hilbert space T-Ld and is thus an 
operator of rank d. Our optimization problem can be 
expressed as follows: 

N 

maximize P c ^J^Tr (pjllj) = 1-Q-P e , (2.2) 
j'=i 

N 

subject to Q = Tr (pH ) = const, n, = /. (2.3) 

3=0 

In Eq. (|2.2[) we introduced the overall error probability 
P e . Since Q is fixed, the maximum of the absolute rate 
of correct results, P c , determines also the maximum of 
the relative rate of correct results at the same value of 
Q, defined as R c = P c /(1 - Q) H|. We thus get 

P™ ax \ Q = R™ ax \ Q (l-Q). (2.4) 

In the special case where R™ ax \Q — 1 errors do not occur 
and the states are unambiguously discriminated. 

As will be shown in Sec. II B, if Q exceeds a certain 
minimum, denoted by Q', the maximum relative rate of 
correct results, R™ ax \Q, stays constant with growing Q. 
The optimum measurement resulting from Eqs. (|2.2[) and 
(I2.3[) with Q > Q' then yields a maximum of P c that is 
smaller than P™ ax \Q> and is therefore without practi- 
cal relevance. One could modify the optimization prob- 
lem and ask for the maximum of P c under the constraint 
that Q does not exceed a certain fixed margin, Qm- The 
maximum of P c under this modified constraint then fol- 
lows from the solution of the original maximization prob- 
lem, Eqs. and and is given by P™ ax \ QM if 
Qm < Q' and by P™ ax \ Q , if Q M > Q' ■ 

In order to derive analytical solutions of the optimiza- 
tion problem posed by Eqs. (|2.2p and (|2.3[) we use the 
operator conditions [25|, [26[ that have to be fulfilled in 
the optimum measurement. Let us begin by re-deriving 
these optimality conditions. For this purpose we intro- 
duce a Hermitian operator Z and a scalar real amplifier 
a. Due to the two constraints in Eq. (|2.3[) the equation 

N 

TrZ-aQ-P c = Ti [(Z - ap)Il } + Tr \{Z - Vjpj)^lj} 

(2.5) 
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is identically fulfilled for any operator Z and any mul- 
tiplier a. Since the detection operators are positive it 
follows that the positivity conditions 

Z - ap > 0, Z - rjjPj > (2.6) 

(j = l,...,N) imply that Tr Z - aQ - P c > 0. Hence 
when the positivity conditions in Eq. (|2.6p are fulfilled, 
the minimum of Tr Z — aQ establishes an upper bound 
for P c . In the optimum measurement, where P c is equal 
to this bound, both sides of Eq. (I2.5[) vanish for the 
optimum multiplier a and the optimum operators Z, ILq 
and Ilj (j — 1,...,N). Due to Eq. (|2.6I) this requires 
that each single term on the right-hand side of Eq. (|2.5[) 
vanishes separately, that is 

(Z-ap)U =0, (Z- Vj p j )U j =0 (2.7) 

(j = 1, . . . , N). Equations jUBJ and ([27f)l together there- 
fore establish sufficient optimality conditions, first de- 
rived by Fiurasek and Jezek [25[ with the help of La- 
grangian multipliers. Using methods of semidefinite pro- 
gramming, the optimality conditions have been shown to 
be not only sufficient, but also necessary (26|. For solving 
the optimization problem it is sometimes advantageous 
to introduce transformed operators [H, [H, [25[ . With 

H = p^Hop 172 , T = p-^Zp^/\ (2.8) 
= P ^^l\ ~ Pj = p-^v 3 p ]P -^ (2.9) 

for j = 1, . . . , N, Eqs. (US]) and §Ff\ yield the optimality 
conditions 

r - Pi > 0, (T - ~ P] )% = (.7 = 1,..., iV), (2.10) 
T-aI>0, (r-a/)n o = 0, (2.11) 

where n + Ylj=i Hj — P> due to the completeness re- 
lation of the detection operators. This representation 
of the optimality conditions results by multiplying Eqs. 
(|2.6j) and (|2.7|) from the left and right by the Hermitian 
operator p -1 / 2 , taking into account that for any opera- 
tor A in Hd the relation Ap^ 1 = can only hold when 
A = 0, since the support of p is the full Hilbert space Hd- 
Provided that an operator Z, a scalar multiplier a 
and positive detection operators IIi,...IIy with n = 
I - X)jL x > satisfy Eqs. {HBJ and (|27H) . or Eqs. (|2~g)) 
- (|2.1ip . respectively, then the detection operators de- 
termine the optimum measurement, which maximizes P c 
with the fixed value Q = Tr (pllo) = Trffo and yields 

JV N 

p c ma 1 Q =E^ Tr (^) = E Tr (pAo- ( 2 - 12 ) 

By taking the trace in both equalities in Eq. (|2.7[) and 
summing over all states in the second equality we arrive 
at the expressions [25| 

Tr (Zn ) =aQ, P™ aa; | Q =Tr[Z(/-n )] =TrZ- aQ. 

(2.13) 



If Q = 0, that is, if inconclusive results are not allowed, 
the optimum measurement corresponds to minimum- 
error discrimination, described by the well-known opti- 
mality conditions arising from Eqs. (|2.6[) and (|2.7[) with 
Ho = and a = 0,0. 

B. Limiting case of sufficiently large Q and relation 
to maximum-confidence discrimination 

When the fixed failure probability Q is getting larger 
and larger, the operator n , defined in Eq. (|2.8|) . will 
turn into an operator of rank d as soon as Q reaches a 
certain minimum value Q' [25| . due to the fact that for 
Q = 1 we must have Hq = I and thus Ho = p with 
rank(p) = d. This means that for Q > Q' the operator 
fto has its support in the full Hilbert space Hd and the 
equality in Eq. (|2.1ip thus can only hold when T—al = 0. 
For Q > Q' the optimality conditions therefore reduce to 

a/-ft>0, {aI-p j )n j =Q (.7 = 1,..., AO (2.14) 

25] . If Q < 1 at least one of the operators ftj has to be 
different from zero. The equality condition in Eq. (|2.14[) 
implies that for any j where lij ^ the eigenstates of 
al—pj belonging to nonzero eigenvalues cannot span the 
full Hilbert space Hd- Therefore when Ilj ^ at least 
one of the eigenvalues of al — fij is equal to zero in the 
optimum measurement with Q > Q' . Together with the 
positivity condition in Eq. (|2.14l) this requires that 

a = max{Ci, . . . Cat}, where Cj — max{eig(pj)}, 

that is, where Cj is the largest eigenvalue of pj. For any 
state j with Cj < a the operator al—fij has its support in 
the full Hilbert space and the equality in Eq. (|2.14p can 
only hold when Bj = 0. For those states j where Ilj ^ 0, 
the support of the operator Ilj is the eigenspace of pj 
belonging to its largest eigenvalue, since this guarantees 
that flj is orthogonal to al—pj j2g. By taking the trace 
in the equality in Eq. (I2.14[) . summing over all states and 
inserting the value of a we arrive at 

Pr x \n = max{Cj}(l -Q) if Q > Q' ■ (2.16) 

In order to determine Q', we have to minimize Q = 1 — 
J2j Tr flj on the conditions that the operators flj have 
the required supports, as described after Eq. (|2.15p . and 
that ffo = p - J2" =1 flj > 0. 

The eigenvalues Cj introduced in Eq. (|2. 15|) have a 
definite meaning in state discrimination. They determine 
the maximum confidence [l6[ that can be achieved for the 
individual measurement outcome j or, equivalently, the 
maximum achievable ratio between all instances where 
the outcome j is correct and all instances where the out- 
come j occurs. In fact, with Eq. (|2.15l) it follows that 

u 3 \ Tr (pUj ) / n, \ TrH, J J ' 1 ' 
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where for each individual j the maximization is per- 
formed with respect to all choices for the detection oper- 
ator Ilj, or Hj, respectively. Clearly, in a measurement 
with Cj = 1 the state j is unambiguously discriminated. 
From Eq. (12.171) it becomes obvious that the maximum 
confidence for the outcome j, equal to Cj, is obtained 
when Ilj has its support in the eigenspace of pj belong- 
ing to the eigenvalue Cj [H, 0, |22| . When this condition 
holds for each of the N states, that is, when for each out- 
come j the confidence is maximal, the measurement is 
called a maximum-confidence measurement [l6j . 

Optimized maximum-confidence discrimination is 
achieved by the specific maximum-confidence measure- 
ment where the probability of inconclusive results takes 
its smallest possible value, Q%g [H, SU. A compari- 
son with the constant Q' characterized after Eq. (|2 . 16[) 
reveals that 



Q' = Q 



MC 
rain 



if Ci = • • • = C N = C. 



(2.18) 



Hence in all cases where the maximum confidence is the 
same for each conclusive outcome, the measurement max- 
imizing P c when Q is fixed at the value Q' is equal to 
the measurement for optimized maximum-confidence dis- 
crimination. When C — 1, the latter measurement cor- 
responds to optimum unambiguous discrimination. 



III. TWO MIXED QUBIT STATES 
A. Method for applying the optimality conditions 

In this paper we want to determine the optimum mea- 
surement for an arbitrary value of the fixed failure prob- 
ability Q, restricting ourselves to cases that allow a sim- 
ple analytical solution. Let us start with the discrimina- 
tion of two mixed qubit states in a joint two-dimensional 
Hilbcrt space, where we use the optimality conditions in 
the form of Eqs. (prTU|) and (|2~TTj) . For N = 2 the trans- 
formed density operators p\ and p 2 defined in Eq. (|2.9|) 
have identical systems of eigenstates, due to the relation 
Pi + P2 = I ■ Their spectral representations therefore can 
be written as 



Pi = Ci|i/i)(i/i| + (1 
h = (l-Ci)h/i)H 



r-C 2 h)H (3.2) 



[18| . where C\ > 1/2 and C 2 > 1/2, since the case 
C\ = C 2 = 1/2 would imply that p\ — p 2 . The con- 
stants C\ and C 2 have the meaning of the maximum 
achievable confidence for the two respective outcomes, 
cf. Eqs. (|2.15[) and (|2.17l) . It is convenient to use the 
orthonormal eigenstates \v\) and \v 2 ) as the basis states 
for solving the optimization problem. For this purpose 
we define the matrix elements 



{vM v i) 
{^\T\ Vj ) 



Pjj 1 



P12 : 
r 12 



\Pl2 

\T12\e 



; (3.3) 

lS (3.4) 



with j = 1,2, where pn + p 22 = 1- In the special case 
when the two given states are pure, pj = \ipj)('ipj\, we 
get from Eq. I|2.9p pj = with p 1/2 \v.j) = y/rjj\il>j), 

which yields 

P12 = VV1V2 (ipi\ip2), Pjj=Vj, Cj = l. (3.5) 

The relation C\ = C 2 = 1 reflects the fact that each of 
the two pure states can be unambiguously discriminated. 
Using Eq. (|2.9[) we obtain Tr (ppj ) = rjj, leading to 

CiPii+Cl-G) p 2 2 = f?i, O--C1) pn+C 2 p22=m- (3-6) 

Now we are prepared to apply the optimality condi- 
tions. We start with the assumption that for the given 
value of Q these conditions are satisfied by a solution 
where Y — al > and where all three detection operators 
are different from zero. Equation (|2.1ip then requires 
that the operators Y — al and IIo have both the rank 
1 and are mutually orthogonal. Likewise we conclude 
from Eq. (|2.10[) that the operators Y — pj and YLj are 
mutually orthogonal rank-one operators for j = 1,2. In 
fact, if for instance IT would be an operator of rank two, 
the equality condition in Eq. (I2.10p could be only sat- 
isfied when r = pi, which would violate the positivity 
constraint Y — p 2 > 0- We are thus led to the ansatz 

Y-aI= IA0XA0I, S = Q|7fo)(ffo|, (AoN) = 0, (3.7) 
Y- pj = Sj = Pj\TVj){nj\, (AjIttj) = (3.8) 

(j = 1,2), where we introduced positive constants (3\ and 
/?2 ■ Here the state vectors characterized by a tilde are 
non-normalized, while 1 7To) is normalized to unity in order 
to yield Trflo = Q. Eqs. d3T7j) and fl3j§ imply that for 
each of the three positive operators Y — al, Y — p\ and 
r — f>2 one of the two eigenvalues is equal to zero. Using 
the orthonormal basis l^)}, this requires that 

|ri 2 | 2 = (r n - a )(r 22 - a) = (r n - Ci)(r 22 -i + c 2 ) 



(r u -i + Ci)(r 22 -c 2 ), 



(3.9) 



due to the fact that the determinants of the three opera- 
tors have to vanish. With the help of Eq. (|3.9[) we obtain 
for each of the three operators the respective eigenstate 
that belongs to its zero eigenvalue. Taking into account 
that these eigenstates are orthogonal to the states \ftj) 
(j = 0, 1, 2) and are therefore proportional to |7fo), |tti) 
and |7r 2 ), respectively, we find that 



TO 



Vr 22 - a\v\) - e 4 Vrn - a \v 2 ) 



(3.10) 



VTTi + r 22 - 2a 

|tti) = y/Tn-l + CM - e- lS ^Y 11 -C 1 \iy 2 ), (3.11) 
|tt2> = V?22 - C 2 Wi) - e^Vrn-l + Cik). (3-12) 

When the expressions under the square-root signs are 
positive and Eq. (|3 .9[) holds, the positivity conditions in 
Eqs. (|2~10l) and ([2~TTj) are satisfied. 

Next we invoke the completeness relation Y^j=o -Ti = 
/. Due to Eq. (I2.8[) the latter relation takes the form 
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n + III + n 2 = p, leading with the help of Eqs. (137)) 
and (|3.8p to the matrix representation 

2 

Pij = Q(Vi\TTo)(*oWj) +^ /3k(Vi\Kk) (KkWj) (3.13) 

fc=l 



for i,j = 1,2. Making use of Eqs. (IXTU|) - (|3~T^1) . we 
arrive at 



Pn + P22 



pu - P22 



\Pi2\e 



i = Q + (/?i+/? 2 )(rii+r 22 -i) 

+ (Px-p2){C 2 -C 1 ), (3.14) 

(rn + r 22 )Q 



rn + r 22 -2a 

+(Pi-fa)(C 1 + C 2 -l) 

jsf Q 



|r 12 |e 4 



+ (/3i + /3 2 )(r 22 -rn) 

(3.15) 

+ /3i + /3 2 Y(3.16) 



T 22 - 2a 



Because of Eq. (|3.9[) we can express Tn, T 22 and |Ti 2 | 
as functions of C\, C 2 and a. Moreover, the positivity of 
/3i and /3 2 requires that e 



as becomes obvious 



from Eq. (ETL6)) . Equations (l37T4j) - (|3~T6]) then repre- 
sent a system of three coupled equations for the unknown 
parameters /3 2 and a. The solution depends on Q. It 
determines the optimum measurement provided that Q 
belongs to a certain region where the expressions under 
the square-root signs in Eqs. (|3. 10[) - (|3.12[) and the re- 
sulting constants /3i and /3 2 are positive. For those values 
of Q where a solution fulfilling these conditions does not 
exist, we have to search for the optimum measurement 
anew. Assuming that one of the two operators Hi or H 2 
vanishes and again supposing that T — al > 0, we obtain 
another solution of the optimization problem in depen- 
dence of Q. In order to determine its range of validity, we 
again have to use the positivity constraints. For the val- 
ues of Q where also this new assumption does not yield 
the optimum measurement we have to drop the supposi- 
tion r — al > 0, that is, we have to consider the limiting 
case r — al = 0. Before proceeding, we separately deal 
with this case. 



B. Solution for sufficiently large Q 

Let us apply the general treatment of the limiting case 
r — al — to the discrimination of two mixed qubit 
states. First we assume that C 2 > C\. As follows from 
the considerations that led to Eq. (|2.16[) . in this case the 
optimality conditions are satisfied when a = C 2 , II i = 
0, and when H 2 = p — ILq is proportional to |^ 2 )(^ 2 |. 
From Trn 2 = 1 - Q we then get ft 2 = (1 - <2)|f 2 )(^2|- 
The constraint IIo = p — II 2 > is fulfilled when the 
determinant resulting from the matrix representation of 
Ho is not negative. This yields the requirement pn(p 22 — 
1 + Q) > |pi 2 | 2 , which holds true for Q > Q\. Here and 
in the following we use the abbreviations 



PU 



IP12T 

Pu 



P22 



\£12£ 
P22 



(3.17) 



Since analogous considerations also apply for C\ > C 2 , 
Eq. (|2.16p can be specified, yielding the maximum prob- 
ability of correct results 



P. 



C 2 (l-Q) for Q>Q' = Qi ifC 2 >Ci, 
'« \Ci(l-Q) for Q>Q' = Q 2 if C 2 <C 1 . 

(3.18) 

When Q = Q' one of the two eigenvalues of IIo vanishes. 
Hence the detection operator Ho is a rank-one operator, 
and since the ranks of Hi and n 2 are and 1, or 1 and 
0, respectively, the measurement is projective. 

In the remaining case C\ = C 2 = C we can again 
apply the considerations that led to Eq. (|2.16p . We then 
conclude that the optimality conditions together with the 
constraint Ho = p — tli — fl 2 > are satisfied when a = C 
and when for j = 1 , 2 



n, 



GjW^iVjl witn (PU - ai)(P22 - "2) > 1/512 1 



(3.19) 

Here the two constants ai and a 2 both have to be non- 
negative, with ol\ + a 2 = Tr (Hi + n 2 ) = 1 — Q. The 
following cases have to be distinguished: 

(i) If |/0x2 1 < niinjpn, p 22 }, the inequality in Eq. (j3. 19[) 
holds true whenever < a.j < pjj — \pi2\ {j — 1, 2), that 
is whenever the given failure probability falls in the range 
2|/?i2 1 < Q < 1- When Q = 2|pi 2 | all three detection op- 
erators have the rank 1. 

(ii) If |pi 2 | > min{pn, p 22 } we first assume that pn < p 22 
and therefore |pi 2 | > p\\. Putting a± = 0, the constraint 
expressed in Eq. (|3 . 1 9[) takes the form we discussed be- 
fore Eq. (|3.17p . which is satisfied if Q > Qi, where the 
measurement is projective for Q = Q\. Similar consider- 
ations hold for p 22 < pn. In summary, for C\ = C 2 = C 
we obtain 




)MC 



< min{pn,p 22 }, 

> Pu, 

> P22- 



(3.20) 
(3.21) 



In Eq. (|3.20p we took into account that for C\ = C 2 the 
value of Q' determines the minimum failure probability 
Qmin necessary for maximum-confidence discrimination, 
see Eq. (|2.18p . We remark that in our earlier paper [l^| 
we calculated Q^fin fo r tw0 mixed qubit states and found 
that the latter is given by the expressions in Eq. ( I3.2ip 
for arbitrary values of C\ and C 2 . 



As outlined in Sec. II after Eq. (|2.4p . a discrimination 
measurement maximizing P c when Q is fixed at a value 
Q > Q' is without practical importance. In the present 
case such a measurement would be a generalized mea- 
surement realizable in an extended Hilbert space, where 
projections onto two orthogonal directions would indicate 
an inconclusive result since Ho has the rank 2, which im- 
plies that also Hg is a rank-two operator. 
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C. Complete solution for the case Ci = C2 = C 

In order to treat the case T — al > we use the method 
developed in Sec. Ill A. It turns out that for C\ ^ C 2 
this leads to a fourth order polynomial equation in the 
variable a(X — a). To obtain a simple analytical solution 
we therefore assume that C\ = C 2 = C, restricting our- 
selves to cases where the maximum achievable confidence 
is equal for the two possible outcomes. Here we present 
the final result for the complete solution, while the de- 
tails of the derivation are shown in Appendix A. Two 
cases have to be considered: 

(i) If I P12 1 < rninjpn, P22} we find that the maximum 
probability of correct results with arbitrary values of the 
fixed failure probability Q is given by 



pi; 



p c (0) (0) 
Q \C(l-Q) 



if < Q < 2\p 12 \, ( ~ 99] 
if 2\p 12 \ < Q < 1, {6 - Zl) 



where 



p(°) = 



1-Q 2C-1 



v / (l-2|p 12 |)(l + 2|p 12 |-2Q). 



2 2 

(3.23) 

(ii) If \pi2\ > min{pn, P22} we assume without lack of 
generality that pn < P22, which means that | P12 1 > Pu- 
Then we get 



Pi 0) (Q) 



if < Q < Q cr , 

if Qcr < Q < Qi, ( 3 - 24 ) 
if Qi < 



,(7(1 -Q) if Qi <Q< 1 
where Qi is defined in Eq. (|3.17|) . Q cr is given by 
2A 



Qc 



with A = P11P22 - IP12 1 (3.25) 



and denotes a critical failure probability that separates 
the regions where the optimum measurement is a gener- 
alized measurement and where it is projective, see below. 
Note that A > since p = r\\p\ + 772/22 is a mixed state. 
In Eq. (I3.24j) we introduced 



P« = C(l-Q) 



2C-1 r 



1-4A 



A(l-2p n ) 



(3.26) 



-(2A-p u )(l-Q) - 2| Pl2 |VA(g-Q 2 -A) 



The corresponding result for p 22 < Pn is obtained when 
in Eqs. (|3.24l) and (|3 . 26[) Q\ and p\\ are replaced by Q2 
and p 22 , respectively. 

When the conditions in the upper lines of Eqs. (|3.22[) 
or (|3.24[) . respectively, are fulfilled, the optimum detec- 
tion operators are given by Eq. (IA8|) together with Eq. 
(fA"4"l) and Eqs. $M\i - (|A11|) . see Appendix A. They de- 
scribe a generalized measurement if Q ^ and Q ^ Q cr - 
For Q = the optimum measurement is equal to the 
projective measurement for minimum-error discrimina- 
tion, yielding the maximum probability of correct results 

p c nax \ Q=0 = p c ME = U 1+Tr \\v2p3 - will) QI- 



1.00 



0.98 



0.96 



<« 0.94 



0.92 



S 0.90 




0.0 
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fixed rate of inconclusive results 
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FIG. 1: Maximum relative rate of correct results 7?™ al |Q = 
pmax^i ^ _ vergus fixed rate of inconclusive results Q 
for discriminating two pure states with 151 = |(^i|'02)| = 0.3, 
0.6, and 0.8, occurring with the prior probabilities 771 = 0.2 
and 772 = 1 — 771 . The dashed lines refer to the regions where 
the optimum measurement is projective, detecting only the 
second state (see text). 



When the middle line of Eq. (|3.24j) applies, the opti- 
mum detection operators follow from Eq. (|A14[) together 
with Eqs. (|A15[) and (j A19|) . In this case the optimum 
measurement is a projective measurement, where only 
one of the states, here the state p 2 , is conclusively dis- 
criminated with a certain probability, while in the pres- 
ence of the other state an inconclusive result is always 
obtained. The lower lines of Eqs. p. 221) or p. 241) . respec- 
tively, correspond to the limiting case Q > Q' discussed 
in the previous paragraph. 

In the special case when C = 1, pn = r)i, P22 = 772, and 
P12 = \/77i77 2 (tpi\ip2) our solution refers to the discrimi- 
nation of two pure states, see Eq. (|3.5|) . Two arbitrary 
qubit states that are mixed can be represented as 



I -Pi 



P2=P2\4>2)('lp2 



1— P2 



2 ' r ' "< T "'^-' 2 

(3.27) 

where and \tp 2 ) are normalized pure states with the 
overlap (ipil^) = S and where Tr7 = 2. The param- 
eters pi and P2 with < pi,P2 < 1 are related to the 
purities of the two states. The condition C\ = C2 = C 

with C < 1 is fulfilled provided that (?7f 1 — l) = rE^it 
see Appendix A. Clearly, for any two mixed qubit states 
pi and P2 there exist particular prior probabilities of oc- 
currence for which our solution holds, and for any given 
prior probability 771 we can find a whole class of state 
pairs where the solution applies. In Figs. 1 and 2 the 
ratio K'l iax \ Q = P c maT | Q /(l - Q) is plotted for the dis- 
crimination of pure states (C = 1) and of mixed states 
(C* < 1). The smallest value of Q where R n c iax \ Q = C 
is equal to the value Q' given by Eq. (|3.21l) and corre- 
sponds to the minimum failure probability necessary for 
unambiguous discrimination if C — 1 [6J , or necessary for 
maximum-confidence discrimination if C < 1 [l8j . 

We emphasize that for Q = and for sufficiently large 
values of Q the analytical solution of our optimization 
problem is known for two arbitrary qubit states, occur- 
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fixed rate of inconclusive results 

FIG. 2: Same as Fig. 1 but including mixed qubit states, 
where C < 1. The parameters are pn = 0.2, |pi2 1 = 0.36, 
and (a) C = 1, (b) C = 0.9, and (c) C = 0.8. When the 
states are represented by Eq. (|3.2Tp . this corresponds to the 
cases that (a) {171 = 0.2, pi = p 2 = 1, \S\ = 0.9}, (b) {qi = 
0.26, pi = 0.92, p 2 = 0.99, l^l = 0.95}, and (c) {771 = 
0.32, pi = 0.90, p 2 = 0.98, l^l = 0.96}. 



ring with arbitrary prior probabilities. The first case cor- 
responds to minimum-error discrimination [l], 0] , where 
for two states the measurement is projective and can be 
easily determined analytically when the Hilbert space is 
two-dimensional. The second case was treated in Sec. 
Ill B and applies when Q > Q' , where Q' is given by 
Eqs. ([3~T5| or ([3~2"Tj) . respectively. For < Q < Q' , how- 
ever, we obtained a simple analytical solution only for 
those states and prior probabilities where the maximum 
confidence is equal for the two outcomes, C\ = C2 = C. 



IV. N SYMMETRIC STATES 

A. Optimality conditions for equiprobable states 
and the limiting case of sufficiently large Q 

Another class of analytically solvable cases refers to 
discriminating N states that are symmetric, which means 
that for j = 1,...,N 



Pj 



= V( j -V piV^-V with V^V = V N = I (4.1) 



[32| . where without lack of generality an arbitrary state 
of the given set of states can be chosen as the reference 
state p\. 

Assuming that each state is prepared with the same 
prior probability, r)j = 1/N, we find that p = 

jfY^j=\Pj — VpV^ and thus [V, p] — 0, which implies 
that the Hermitian operator p and the unitary operator 
V can be diagonalized in the same orthonormal basis [32| . 
Introducing the orthonormal basis states |n), we get the 
spectral representations 

N d d 

vEft=E r 'l r ^ f 'l' V = J2vi\n)(ri\, (4.2) 



N 



where = NS W and \ Vl \ 2 



,.N 



1, see also 



[22l | . We now focus on the optimum measurement, where 
the optimality conditions, given by Eqs. (|2.6p and (|2.7|) . 
are satisfied. Let us suppose that ill is an element of the 
set of optimum detection operators. We introduce 



(4.3) 



U = 1, 



. , N), where due to the symmetry the operator 
_ ( _ ,~j = I — n commutes with V. Using [V, p] = 0, 
this implies that 



y N n- 



[n , v] = po,p] = [n , z) = [ P , z\ = [v, z] = o. (4.4) 

Here the second equality sign follows from the fact that in 
the optimum measurement (Z — ap)Ho = Hq(Z — ap) = 0, 
which leads to [IIo,.Z] = a [Ho,/?]. Upon inserting Eq. 
(|P|) into Eqs. (I^d| and (|2~71) with rjj = 1/N, taking 
Eq. (|4.4I) into account, it becomes obvious that if ill ful- 
fills the optimality conditions, then these conditions are 
also fulfilled by each of the operators IL; defined in Eq. 
(|4.3p . in analogy to the derivation of the optimality con- 
ditions for maximum-confidence discrimination of sym- 
metric mixed states in our previous paper 22[. Hence 



the detection operators for the optimum measurement 
can always be chosen in the form of Eq. (|4.3|) . The opti- 
mality conditions thus reduce to the conditions 



Z 

z- 



ap > 0, (Z -ap) n =0, (4.5) 



which were first derived by applying group-theoretical 
methods [26|. Using Eq. (14.31) and the properties of V 
we find that 



n n 



d 

E 

i=i 



d 

E 

i=i 



zi\n)(ri\ 



(4.7) 

where we introduced the spectral representation of Z, 
taking into account that n and Z commute when the 
optimality conditions are fulfilled. With the help of Eqs. 
(|2.12p and (|4.3[) we can represent the maximum proba- 
bility of correct results achievable with the fixed value 
Q = Tr(pn ) as 



P" 



\Q 



Tr(pini) where Q = 1 - APTr (pIIi). (4.8) 



Let us discuss the limiting case of sufficiently large Q, 
which we treated in general in Sec. II B. According to 
Eq. (|2.14l) in this limit the optimality conditions, Eqs. 
(|4.5[) and (|4.6[) . reduce to the conditions al — p~\ > and 
(al — pi) LTi = 0, where we used the transformed oper- 
ators introduced in Eqs. (|2.8p and (|2.9p . The spectral 
representation of p~\ can be written as 



Pi 



1 



-1/2 



PIP 



-1/2 



3 = 1 



1 = 1 



k\ d 

cEk>H+ E v M){»k\ 

k=l k=ki + l 

(4.9) 
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where C denotes the largest eigenvalue and ki is its de- 
gree of degeneracy. The optimality conditions are fulfilled 
when a = C and when the support of Hi = p^^Iiip 1 / 2 is 
the Hilbert space spanned by the eigenstates \ v\) . . . \ vki)- 



We then obtain P" 



\Q 



CNTrlli, which yields 



= C{1 -Q) if Q>Q' = Q 



MC 
rain ' 



(4.10) 



Here we took into account that the maximum confidence 
is the same for discriminating each of the equiprobable 
symmetric states pj = V^^ 1 ' piV^- 1 ) since the largest 
eigenvalues of the operators pj are identical, that is, 



C\ — • • • — Cjv — C. 



(4.11) 



Q' therefore corresponds to the failure probability Q^E 
in optimized maximum-confidence discrimination, or, if 
C = 1, in optimum unambiguous discrimination, re- 
spectively, see Eq. (|2.18[) . In order to determine Q' 
we have to find the minimum of Q = 1 — TVTrlli 
on the condition that with the given support of III 
the positivity constraint Ho > is satisfied, which 
because of Eq. (|4.7p is equivalent to the constraint 
Ya=i [ r i ~ N {r^IV^n)] \n){n\ > 0, see [Hj]. In general, 
this is a nontrivial task. In the special case when the 
largest eigenvalue of p\ is nondegenerate, TLi is propor- 
tional to | i^i ) (i/i | . In agreement with our previous paper 

p2j we then obtain Q' 
ki = 1. 



' rain 



1 



t{k^>f} 



for 



B. TV equiprobable symmetric pure qudit states 

1. Solution for a special class of states 

First we consider the discrimination of N symmet- 
ric equiprobable pure states spanning a ci-dimcnsional 
Hilbert space T-Ld- This means that d < N and implies 
that all expansion coefficients q of the states with respect 
to the d-dimensional eigenbasis of the symmetry operator 
V are different from zero (22[. In general, N symmetric 
pure qudit states are given by 



|Vi) = Vtt-%i) with 



d 

i=i 



Cl \n) (j = l,...,N). 

In order to obtain an analytical solution we restrict our- 
selves to special states where only two different values of 
the expansion coefficients q occur. We assume that 



hM=ci^|n)+C2^ In), m\ Cl \ 2 + (d-m)\c 2 \ 2 = l, 

1=1 l=m+l 

where m < d/2 and where we distinguish the cases 
|ci| > | C2 1 and |ci| < \c 2 \ with ci,C2 ^ 0. In Appendix 
B we derive the maximum probability of correct results 



P™ ax \ n with the fixed probability Q of inconclusive re- 
suits, making use of the optimality conditions, Eqs. (|4.5[) 
and g^]). Using Eq. (|BT0|) with p = 1 and Eq. (|BT3l) 
we obtain 



PI: 



with 



P(°): 



P { c\Q) if Q < Q' =l-dmin{|ci| 2 , 
£(1-Q) if Q>Q' 



\C2?} 

(4.14) 



myjm\ci\ 2 - Q + (d- m)\c 2 \ 



|c 2 | 



>1 



'ifH<i. 

|c 2 | - 



Vd-my/(d-m)\c 2 \ 2 - Q + m\ci\ 

' (4.15) 

The lower line of Eq. (|4.14l) follows from Eq. (14.10j) and 

from the fact that P C (0) (Q')/(1 - Q') = C = 4. The 
optimum detection operators are given by Eqs. (|B3[) and 

(|B71). 

When Q = the optimum measurement corresponds 
to minimum-error discrimination. For this case the op- 
timum detection operators discriminating the equiprob- 
able states given by Eq. (|4.12[) are known to be Hj IE = 

jjP~ 1 ^ 2 \i'j)(ipj\p~ 1 / 2 [Hf, yielding for an arbitrary ref- 
erence state |i/>i) the maximum probability of correct re- 
sults Pc nax \ Q=0 = w fef=i Nl) • On the other hand, 

when Q — Q 1 — Qmi n the optimum measurement cor- 
responds to optimized maximum-confidence discrimina- 
tion. For N equiprobable symmetric pure qudit states 
with an arbitrary reference state it has been shown 
that Q™g = 1 - dmim{|c/| 2 } [H|,[22j], including the spe- 
cial case of unambiguous discrimination when N = d 
[33j . The results obtained from Eq. (|4.14l) for the max- 
imum confidence C, for Q' — Qjf?, and for P™ ax \q =0 
agree with these previous results. Note that for |ci| = 
(cal = 4 we obtain Q' = 0, which means that in this case 
minimum-error discrimination and optimized maximum- 
confidence discrimination are equivalent, cf. also Sec. V 
A. 

We emphasize that while in the limiting cases Q = 
and Q — Q' the optimum measurement with a fixed value 
of Q is known for N equiprobable symmetric pure qudit 
states which are arbitrary, our complete analytical solu- 
tion, Eq. (14.14[) . that interpolates between these limiting 
cases, is restricted to the special class of those symmetric 
states where the reference state is given by Eq. (|4.13l) . 



2. Application to N linearly independent symmetric states 

When the number of pure states is equal to the dimen- 
sion of the Hilbert space spanned by them, N — d, the 
states are linearly independent. This means that without 
lack of generality the eigenvalues of the symmetry oper- 
ator V in Eq. (|4.2|) can be written as vi = exp (2iri l -j^ 
for I = l,...,iV~j33|. With the help of Eqs. (|4TT2"j) 
and (|4.13l) we obtain for d = N the mutual overlaps 
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(V>j#fe) = Ei=i v i v * 3 \ l CPA r l)\ 2 ^ yielding with j ^ k 

m 

(^# fe ) = (|ci| 2 -| C2 | 2 )^exp[27r^(fc-i)]. (4.16) 
l=i 

Clearly, when m — 1 all mutual overlaps are equal and 
real. If d = N = 3 this is always the case, because of 
the requirement m < d/2 in Eq. (|4.13|) . For d = N > 4 
and m > 2, however, different values of the overlap occur 
for the linearly independent symmetric states specified 
byEq. KWi . 

We now specialize to N linearly independent symmet- 
ric pure states with equal mutual overlaps S, where 



(^ k ) = S = \ Cl f 



C2 



= 1-N\c 2 \ 2 



(4.17) 



according to Eq. (|4.16l) with m = 1. Using the normal- 
ization condition given in Eq. (I4.13[) . we get the expres- 



■is 



and 



l+W 

sions \ci\- = — y — 
conclude that — (JV 

|c 2 ' 9 ■ 
into Eq. f|4. 1 5[) we obtain for \c\ 



■1)- 



1 < S < 1 since neither 



from which we 
| C! | 2 nor 

2 can be negative. Upon inserting these expressions 



> \c 2 \, that is for S > 0, 



p(0) = 



1 



JV 

if < S < 1 



1 + (7V-1)5 
JV 

for 



Q + (N 

Q<Q' 



l-S 

N 



s. 



(4.18) 



while for negative overlaps we arrive at 



p(0) : 



1 

JV 



if - 1 v L t<S'<0, for Q < Q' = (N-1)\S\. (4.19) 



If Q > Q' Eq. gHI yields P c maa; |Q = 1 - Q, which cor- 
responds to unambiguous discrimination since the proba- 
bility of errors vanishes in this case, P e = 1 — Q — P c = 0. 
The minimum failure probability necessary for unam- 
biguous discrimination is thus given by Q' . When 5 = 
Eqs. (|4.18|) and (|4.19j) apply only for Q — and we 
get Pc nax \q- = 1) as expected for mutually orthogonal 
states. On the other hand, in the limit where S ap- 
proaches — N _ 1 and Ci therefore approaches zero, the 
states get linearly dependent and span a Hilbert space of 
dimension N — 1 since m = 1. For example, if JV = 3 
and (tpj\i^k) = —\ = cos 2 ^ (j ^ fc) we arrive at the 
trine states that can be represented by three real sym- 
metric state vectors spanning a two dimensional Hilbert 
space. We note that for N = 2 Eqs. (|47L8|) and (|4TT9|) 
are identical and reproduce the result obtained for the 
discrimination of two equiprobable pure states [23j . For 
N > 3, however, the explicit expressions for p™ ax and 
Q' depend on the sign of the overlap between the states, 
see Fig. 3. 
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FIG. 

-ymax 



Maximum relative rate of correct results R™ ax \cj = 
P™ ax \Q/ (1 — Q) versus the fixed rate of inconclusive results Q 
for discriminating N equiprobable linearly independent sym- 
metric pure states with equal mutual overlaps S = —0.3 (full 
lines) and S = 0.3 (dashed lines). For N — 2 both lines 
coincide. 
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FIG. 4: Same as Fig. 3 for discriminating N = 3 equiprobable 
symmetric mixed qubit states, characterized by Eq. (|4.20[) 



with |ciT = 0.85 and 
purity parameter p. 



'■2 



= 0.15, for different values of the 



state \tp) and with < p < 1. Hence the most general 
representation for N symmetric mixed states in a two- 
dimensional joint Hilbert with the identity operator / 
is given by the density operators pj — yW-iJ^ytu'- 1 ) 
with j = 1, . . . , N and with V'V — V N — /, where 



Pi =p\ipi){^i 



l-p 



I, hh>=ci|ri)+C2|r 2 >. (4.20) 



Without lack of generality we assume that |ci| > \c2\- 
Using Eqs. (|B8|) . (fB10|) and (IBM)) (see Appendix B) we 
find that the maximum probability of correct results with 
the fixed failure probability Q is given by 



{Pc 0) (Q) 
\C(1-Q) 



if Q<Q'= P (l-2|c 2 | 2 ) 
if Q>Q' 



(4.21) 



C. iV equiprobable symmetric mixed qubit states 



The density operator of any mixed qubit state can be 
written as p \ip) + I with a certain normalized pure 



with P (0) = Ji 
N 




2Q 



l+p(l-2|c 2 | 2 ) 

Q 



l+p(l-2ca 



K"2 



(4.22) 
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The ratio P™ ax \Q/(l — Q) is plotted in Fig. 4. By calcu- the normalized eigenstates as \ip k ) = ^ z uu\Sj), taking 
lating C — Pc°\Q / )/(1 — Q') we arrive at the maximum into account that the expansion coefficients u a = (i(j k \s lx 



confidence C 



1 



2p Icicaj 



and the eigenvalues Xk are identical for the different 



x /l-p2(l-2|c 2 |2) 

with the result derived for C in our previous paper [22j | 
Q' coincides with the smallest failure probability neces- 
sary to achieve maximum-confidence discrimination [22| 
On the other hand, for Q = we get the result 



, m accordance states j since the matrix elements (s l Apj\s l j) do not de- 



pend on j. Hence for each index k also the N eigenstates 
■ ■ ■ , \ip%), belonging each to a different density op- 
erator, represent a set of linearly independent symmetric 
states. 

Let us introduce the iV-dimensional subspace H k 



P n 



2=o 



P 



ME . 



1 . 

— (l+p|cic 2 | 



(4.23) 



which means that using our general solution for p™ oa: 

we determined the maximum probability of correct re- 
sults in minimum-error discrimination of the N mixed 
qubit states. 



N 

spanned by the eigenstates \ip k ), ■ ■ . , \^%) and let 1^ 
be the identity operator in this subspace. Since the rela- 
tion (if>i\ipj ) — (k 7^ k') holds true for any j, j, the D 

L, . . 



The special case \c\\ 



C2 



0.5 is worth men- 



tioning, where Q' = and therefore P c maa: /(l - Q) = 
C = for any value of Q. The measurements 

for minimum-error discrimination and for optimized 
maximum-confidence discrimination are then the same, 
cf. Sec. V A. When N — 3 this applies for the de- 
polarized trine states, described by Eq. (|4.20[) with 
|^) = ^[exp(if) | ri )+expHf)|r 2 )]. 



D. N special symmetric mixed states of rank D 
spanning a joint Hilbert space of dimension ND 

In the following we treat a special case where the 
optimum measurement for discriminating N symmetric 
mixed states with a fixed probability Q of inconclusive 
results can be obtained by applying the pure-state solu- 
tion. We consider N mixed states of rank D, occurring 
with the prior probabilities rjj (j = 1, . . . , N) and being 
described by the special density operators 



subspaces ~H k N with k — 1, . . . ,D are mutually orthogo- 
nal, which implies that ^feLi 1^ = I n the following 
we use the convention that operators labeled by the su- 
perscript (k) act in the subspace %%■ Introducing the 
symmetry operator with \*fj k ) = V {k ^_ x ), it fol- 
lows that pj = Vpj-iV\ where V = <g) . . . ® V (D \ 
Clearly, the N mixed states given by Eq. (|4.24[) or, equiv- 
alcntly, by Eq. (|4.26[) are symmetric, but in the following 
we do not necessarily assume that they occur with equal 
prior probabilities rjj. 

The optimum measurement, discriminating the mixed 
states with a maximum probability of correct results 
for a fixed value of the failure probability Q, is deter- 
mined by the optimality conditions, Eqs. (|2.6|) and (|2.7p . 
where because of Eq. (I4.26[) p — J2k=i ^ k w ^ n 



P 



(*0 



J2f=i Vj pf- Using the orthogonality of the dif- 
ferent subspaces H% it follows that the optimality condi- 
tions are satisfied for a certain value of the real multiplier 
a and for certain operators Z and II; (i = 0, 1, . . . ,N) 
when 



D 



N 



z = ^x k z^ , il = £ nf 5 with £ nf > = /<*> , 



k=l 



k=l 



D 

ft=E s «'i4x4'i 

2 = 1 



with (s l M}=S l 



(4.24) 



where we assume that the overlaps of basis states be- 
longing to different density operators obey the special 
relation 



(s\\ s l ')=SS w 



with - < S < 1 



(4.27) 

provided that Z^> and Ilf ) {i = 0, 1, . . . , N) satisfy the 
corresponding optimality conditions in their respective 
subspaces 71% for the same value of a. The latter re- 
quirement can indeed be fulfilled in our problem, since 
in each of the different subspaces T-l% the structure of 
the resulting optimality conditions is identical, due to 
the fact that (ipi\ipj) does not depend on k. In par- 
ti ^ j) (4.25) ticular, this means that in the optimum measurement 



Tr^nf 5 ) = Tr(^ ij lt- 1J ) for k = 2,...,D. Con- 
sequently, using again the orthogonality of the differ- 
ent subspaces %%, we find that the detection opera- 
tors given in Eq. (|4.27[) maximize P c at the fixed fail- 
ure probability Q = Tr (pllo) that can be written as 
Q = EfeAfeTr^WlT^) = TrGoWlI^). According to 
Eq. (|2.12p they yield the maximum probability of cor- 
rect results 



and with 1,1' = 1, . . . , D. For each value of I the states 
\s[), . . . , |sjy) form a set of N linearly independent sym- 
metric states with equal mutual overlaps, as becomes ob- 
vious from Eq. (|4.17|) . The spectral representations of 
the density operators can be written as 



Pi=f> fcP f with pf=I^X^I, 



fe=l 



D 

E 

fe=i 



A fc = l, 
(4.26) 

where {^ k \ij} k ')=5 kk > and (ip k \ip k ' ) = S 5 kk > for i^j. The 
latter equation follows from Eq. (|4.25p after expanding 



3 = 1 



k=l 



N 

=E 

3 = 1 



r]jTr {^pf 



(4.28) 
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The dependence of P™ ax on Q is thus exactly the 
same as the respective dependence that arises from 
the corresponding discrimination problem for the states 
),..., \iPn) i n an y one 01 the orthogonal subspaces 
W.% (k = 1, ...,D). Without lack of generality in Eq. 
(|4.28[) we referred to the subspace "H^. In other words, 
provided that the pure-state optimization problem can 
be solved for the set of linearly independent symmetric 
states occurring with the prior probabilities r\j 

(j = 1 . . . , N), that is when the optimum detection op- 
erators in the subspace can be determined, we know 
the complete solution. 

We emphasize that P™ ax \g does not depend on the 
rank D of the mixed states pj nor on their matrix ele- 
ments su>, or their eigenvalues Afc, respectively. Using 
the pure-state results derived in this paper, we obtain 
analytical solutions in two cases: 

(i) For N mixed states that are defined by Eq. (I4.24j) or, 
equivalently, by (|4.26l) and occur with equal prior proba- 
bilities 7/ j — 1/N, the final result P" lax \g is represented 
by Eq. (|4~T4]) with d = N, together with Eqs. (|4~Tg|) and 
(|4"H?|) . The failure probability Q' , given by Q' = S if 
S > and Q' = (N - l)\S\ if < S < 0, now corre- 
sponds to the minimum failure probability necessary for 
unambiguously discriminating the N mixed states. 

(ii) When N — 2 and rji < r) 2 the final result P™ ax \ n 
for optimally discriminating the two mixed states de- 
scribed by Eqs. (|4.24l) or (|4.26|) . respectively, is ob- 
tained by substituting the expressions pn = rji, P22 = V2, 
P12 = y/mmS and C = 1 into Eqs. (|X2^|) - (pOo]) . thus 
using the corresponding solution for the optimum dis- 
crimination of two pure states occurring with arbitrary 
prior probabilities. We note that two pure states can be 
always written as symmetric states with respect to a suit- 
able basis. Our general result for PJ 71031 1 reveals that the 
minimum failure probability Q' required for unambigu- 
ously discriminating the two mixed states is given by Eq. 
(|3.21j) . in agreement with earlier results for the optimum 
unambiguous discrimination of these two special mixed 
states [lll[I2. 



FURTHER APPLICATIONS 



iV equiprobable mixed qudit states resolving 
the identity operator 



When the dimension d of the joint Hilbert space is 
larger than two and the states are genuinely mixed, that 
is when the discrimination problem cannot be reduced to 
the problem of discriminating pure states, it is in general 
hard to obtain analytical solutions. However, there is an 
exceptional case. We consider N states that occur with 
equal prior probabilities, rjj — 1/N , and are described by 



the special density operators pj (j — 1, . . . , N), where 



1-p 1 N 

pj = v 1 + — — 1 with p = — 



^ = ~d- 

3=1 

(5.1) 

Here / is the identity operator in T-ld, and < p < 1. 
Eq. (|5.ip means that the identity operator can be re- 
solved as a weighted sum over the density operators pj. 
The largest eigenvalue of any one of the operators pj = 
p-i/2^p^p-i/2 _ JLpj^ determining the maximum confi- 
dence Cj of the outcome j, is given by (pd+1— p)/N, and 
the corresponding eigenstate is l^,-). Taking into account 
that p is proportional to /, it follows from the consider- 
ations in Sc. II B that for maximum-confidence discrim- 
ination both the operators Ilj and H,- = p~ 1 / 2 fljp~ 1 / 2 
are proportional to \ipj)(ipj\- Since Eq. (|5.1j) implies that 



7T E,=i l^jX^I = a, the operators 



H =0, Uj = -MAW (j = l,...,N) (5.2) 

fulfill the completeness relation. Therefore maximum- 
confidence discrimination is possible without inconclusive 
results, and we get with the help of Eq. (|2.18|) 

Ci = • • • = Cn = C = ^--z -, Q' = Qmii = 0- 



N 



From Eq. (|2.16p we then obtain 



P" 



= C(l -Q) if Q > 0, 



(5.3) 



(5.4) 



which in particular means that P™ a: 
The detection operators given by Eq. 



_ l+p(d-l) 



2=0 



N 

(I5.2[) describe op- 
timized maximum-confidence discrimination since they 
yield the smallest possible failure probability, Q = 0. On 
the other hand, they also describe minimum-error dis- 
crimination, since they maximize the probability of cor- 
rect results when Q is fixed at the value Q = 0. Hence 
both measurements coincide and the relative rate of cor- 
rect results cannot be increased by admitting inconclu- 
sive results. 

Eqs. (|5.2j) - (|5.4[) agree with our previous results [22[ 
obtained by studying the optimized maximum-confidence 
discrimination of symmetric states obeying Eq. (|5.1|) . 
However, they show that these results also hold in a more 
general case. While the derivation in [22j supposes sym- 
metric states \ipj) (j — 1,...,N) with equal expansion 
coefficients q with respect to the basis of the symmetry 
operator, q = 1/Vd for I = l,...,d (see also Appendix 
B), we emphasize that in Eq. (|5.ip the N states pj, or 
\ijjj), respectively, need not necessarily be symmetric. 

Moreover, the expression for P™ ax \q_ Q following from 
Eqs. (15.31) and (I5.4[) generalizes the known result for 
minimum-error discrimination of pure states resolving 
the identity operator [U, [35[ to a special class of mixed 
states. We mention that the corresponding pure-state 
measurement has been experimentally realized for a set 
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of qubit states (d = 2) with N = 3, given by the sym- 
metric trine states, and also for a set with N = 4, given 
by the tetrad states [34] . The latter are defined as 

m) = S= = 1,2,3), m)=\o). 

(5.5) 

Note that the tetrad states do not belong to the class of 
fully symmetric states considered in the previous section, 
but possess only a partial symmetry. It is easy to check 
that they fulfill the requirement \ X)t=i iV'iKV'.jl = f- 
When p < 1 the end points of the Bloch vectors be- 
longing to the four depolarized tetrad states resulting 
from Eq. (|5.1[) together with Eq. (|5.5|) form a regular 
tetrahedron within the Bloch sphere, with circum radius 
p and edge length p-v/8/3. In the special case N = 4 
and e? = 2 we get P™ ax \Q =0 = (1 + p)/4, in accordance 
with the result following from a recent solution [361 ] for 
the discrimination of mixed qubit states with Bloch vec- 
tors forming a regular polyhedron, where an approach to 
study minimum-error discrimination was applied that is 
based on a geometrical method usin g H elstrom families 
of ensembles in convex optimization [36l |37{ . 

B. Partially symmetric states 

The methods developed in this paper can be extended 
to the problem of discriminating between states that pos- 
sess a certain partial symmetry. We assume that the 
given set of states consists of two sets of a equiproba- 
ble symmetric states and that the total density opera- 
tors resulting from each set commute. More precisely, 
we refer to the discrimination of N states (N > M > 1), 
where the states p\ , . . . pm are symmetric as described in 
Sec. IV, occurring with equal prior probabilities rjx/M, 
and where the remaining states pM+it ■ ■ ■ Pn with equal 
prior probabilities (1 — r)i)/(N — M) are also symmet- 
ric. In addition, we assume that the operators VJ . =1 Pi 

and X^Lm+i Pi nave the same eigenbasis, which implies 
that the symmetry operators V and U , referring to the 
two symmetric sets, are both diagonal in this eigenbasis. 
In analogy to the derivation of Eqs. (|4.5I) and (|4.6I) it 
follows that in the optimum measurement the detection 
operators can be supposed to obey the same symmetry 
as the density operators, that is, H, = V^^'TliV^^ -1 ' 
for j = 1,...,M and U M +j = U^-^IL M +iU^^ for 
j = 1,...,N-M. Eqs. (J!!]) and (gTTj then reduce to 
the optimality conditions 

Z - ap > 0, (Z- ap) n = 0, (5.6) 

Z _V2PM ±1> ( z _mPM ±1 \ 

N-M ~ ' V N-M J + v ' 

where 771+772 = !■ These conditions contain only n and 
the two detection operators belonging to the reference 



states p\ and pM+i of the two symmetric sets. 



VI. RELATION TO STATE DISCRIMINATION 
WITH A FIXED ERROR PROBABILITY 

A. General considerations 

So far we considered the measurement strategy that 
maximizes the overall probability of getting a correct re- 
sult, P c , with a fixed value of the failure probability Q. 
The general relation P c + P e + Q = 1 implies that 

Pr X \ Q = l-Q-Pe(Q), (6.1) 

where P e (Q) is the minimum overall error probability 
that can be obtained at the same fixed value of Q. An- 
other discrimination strategy maximizes P c under the 
constraint that the overall error probablity P e has a fixed 
value [28Ll29l|. We then get 

Pr x \ Pc = l-Pe-Q(Pe), (6.2) 

where Q(P e ) is the minimum failure probability necessary 
to achieve the same fixed error rate P e . 

Let us investigate the relation between the optimiza- 
tion problems posed by these two strategies. For this 
purpose we suppose that P e (Q), introduced in Eq. (|6.1|) . 
is a monotoneously decreasing function of Q in a certain 
interval around a value Q — Q a . Then from the assump- 
tion P e (Q a ) = P e Q (1) it follows that Q(P") = Q a (2), 
where Q(P e ) is introduced in Eq. (I6.2[) . In order to verify 
this intuitive statement, we use an indirect proof. Sup- 
pose that Q(P") — Qf) where Qp > Q a . This means that 
the failure probability Q a is not large enough to achieve 
the value P", or, in other words, when Q is fixed at Q a we 
get a minimum overall error probability P e that is larger 
than P e Q , in contradiction to the assumption (1). Now 
suppose that Q(P e Q ) — Qp where Qp < Q a - This means 
that the value P" can be already reached at a value of 
Q that is smaller than Q a which together with the as- 
sumption (1) is a contradiction to the fact that P e {Q) is 
monotonously decreasing. Hence the conclusion (2) in- 
deed follows from the assumption (1). In an analogous 
way the equation P e {Q a ) = P" can be shown to follow 
from Q(P") = Q a . These findings are summarized as 

P e (Q a ) = P e Q o Q(P") = Q a . (6.3) 

Taking Eqs. (|6.1[) and (|6.2I) into account, we thus ob- 
tain the relation Pf" 1 1 _ = p rnax I . Hence we con- 

elude that the detection operators maximizing P c at the 
fixed failure probability Q a (or minimizing P e at this 
value Q a , respectively), are the same as the detection 
operators maximizing P c at the fixed error probability 
= 1 — Q a — p™ ax I ( or minimizing Q at this value 

P" , respectively). 
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The latter conclusion can be also obtained in more 
formal terms. Let us assume that the detection opera- 
tors Ilo, III, ■ • ■ , Iljv, the operator Z and the scalar mul- 
tiplier a fulfill the optimality conditions, Eqs. (|2.6[) and 
(|2.7p . Then the detection operators determine the mea- 
surement that maximizes P c at the fixed failure proba- 
bility Q — a _1 Tr(ZIIo), as becomes obvious from Eq. 
(I2.13[) . This yields the probability of correct results 
pmax\ = i _ p e _ q = t t z - a Q, see Eq. (l2~13l) . 
From the latter equality it follows that in the optimum 
measurement P e = 1 — Tr Z — (1 — a)Q. Taking again 
Eq. (|2.13l) into account, we conclude that the same op- 
timum detection operators also characterize a measure- 
ment which maximizes P c at the fixed overall error prob- 
ability P e = 1 - Tr Z - {a- 1 - l)Tr (ZU Q ). 

Due to the connection between the optimization prob- 
lems in the two strategies we can directly determine the 
solution P™ ax \ Pc from the solution P c nax \Q and vice 
versa, using Eqs. (|6.ip and (16.21) and taking into ac- 
count that according to Eq. (|6.3p the function Q(P e ) 
is the inverse of the function P e (Q). We note that the 



determination of PI' 



is only of practical interest if 



P < PR 



1 - P" 



l<2=o 



(6.4) 



with Pe denoting the minimum error probability ob- 
tainable in the strategy of minimum-error discrimination 
[3; HI > where inconclusive results do not occur. Admitting 
a larger value of P e does not yield any advantage since 
for P e > P E we get P™ ax \ Pc < 1 - P e < 1 - Pe, where 
Eq. (|6.2p has been used. 



B. Example 

In the following we present an example where we de- 
rive the maximum probability of correct results with a 



fixed error rate, P? 



, with the help of the result for 



P™ ax \q- We consider the discrimination of two mixed 
qubit states with C\ = C 2 — C, see Sec. Ill C, assuming 
that pn < p 22 . From Eq. (pT|) together with Eqs. (|3~2"2]) 
and (j?72"4")) we find that 



Pe(Q) 



1-Q 2C - 1 



V(l-2M)(l + 2|pi S 



20), 
(6.5) 

provided that Q is restricted to a certain interval. Insert- 
ing the boundaries of this interval into Eq. (|6.5p leads 
to the restriction P' e < P e {Q) < Pe, where Pe refers to 
minimum-error discrimination, see Eq. (|6.4p . and 



P e (2\ Pl2 \) = (1 - C)(l - 2|pi 2 |) if |pi 2 | < P11, 

Pe(Qcr) if |/3l2 j > Pll 

(6.6) 

with P e (Q cr ) = lB^0fl + (1 - C)(p 22 - Pll ). From 



the function P e (Q) we can derive the inverse function 



<5(P e ). Making use of Eq. (|6.2I) we obtain 

Pr x \ Pc = Pe + (2C-1) 2 (1 - 2\ Pl2 \) (6.7) 

+2(2C-l)V(i - 2|p 12 |)[P e -C(l-C)(l - %!■ 
if P' e <Pe< Pe, 



which reduces to P™ ax \ Pc = (v^+V 1 2 |p 12 |) 2 when 
C = 1, that is when the states are pure, in accordance 
with (29[ . Next we use Eq. j6Tj together with Eq. (|3T20|) 
and get P e (Q) = (1 - C)( l - Q) for Q' < Q < 1, where 
Q' is given by Eq. ([3T5T]) . This restricts P e (Q) to the 
interval < P e (Q) < P e " with 



p': 



{\-C){\-2\p X2 \)=P' e if |pi 2 |<pu, , fi8 ^ 
(1-C)(l-Qi) if |pr,l > - lb ' 8j 



Pn- 



For 1 we can determine the inverse function Q(P e ) 
and insert it into Eq. (|6.2p . arriving at 



= P e 



C 



l-C 



if 0<Pe<P" (C^l). 



(6.9) 

Equation (|6.9p reflects the fact that unambiguous dis- 
crimination, where P e = 0, is impossible when C =/= 1, 
since then p™ ax = which means that Q = 1 and the 
measurement always fails. 



If IP12I < P11 Eqs. 



(|6.9p determine the complete 



solution since in this case P' e = P". For C = 1 the 
solution is in agreement with the result obtained for pure 
states b y d irectly performing the optimization when P e 
is fixed [29j |. 

If \pi 2 \ > P11, that is if Pg 7^ P", we still have to 
consider the interval P e " = P e (Qi) < P e < P e (Q cr ) = P P '. 
From Eqs. and ([AlO)) we find that 

P e (Q) = (1 - C)(l - Q) + (2C - l) 7 i (Q) (6.10) 



if Qcr < Q < Qi, where 71 (Qi) = 0. Us ing E q. (|AT9l 
it is in principle possible to invert Eq. (|6.10p . that is, 
to determine the function Q(P e ), which yields P™ ax \p a 
according to Eq. (|6.2p . Since for C ^ 1 the calculations 
are rather involved, we specialize to the case C = 1, 
where P'J = and P e (Q) = 71 (Q) , which due to Eq. ([672]) 
yields P c maa: |p c = 1 - 71 - Q(7i). The right-hand side of 
the latter equation is equivalent to 72 , as becomes obvious 
from Eq. (j A19|) . Substituting 71 = P e into the equation 
for 72 given by Eq. (j A19|) . we obtain for \pn\ > pn 



if G = 1 and < P e < 



3.11) 



(£1 



-\p12\y _ p i 

-2|pi 2 | r e> 



where A = pnp 22 —\pi 2 \ 2 . Equation (|6.1ip together with 
Eq. (|6.7I) for C = 1 determines the complete solution 
when the states are pure with \pi 2 \ > pn. Using Eq. 
(|3.5|) we find that this solution is in agreement with the 
results obtained for two pure states by directly perform- 
ing the optimization for a fixed value of P e [29J . 
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VII. SUMMARY AND CONCLUSIONS 

In the main part of the paper we considered a measure- 
ment for state discrimination that minimizes the error 
probability P e , or maximizes the probability P c of cor- 
rect results, respectively, when a certain fixed probability 
Q of inconclusive results is admitted. For a number of 
problems not treated before we derived analytical solu- 
tions for the optimum measurement: 

(i) We investigated the discrimination of two arbitrary 
mixed qubit states that occur with arbitrary prior prob- 
abilities. For the case that the two conclusive outcomes 
can be discriminated with the same maximum confidence 
we obtained the complete solution, see Eqs. p.22[) - 
(|3.26p . This solution includes the discrimination of two 
pure states occurring with arbitrary prior probabilities. 

(ii) We studied the discrimination of N symmetric states 
spanning a d-dimensional Hilbert space. For a cer- 
tain class of symmetric equiprobable pure qudit states 
(d < N) we derived the solution, given by Eqs. (I4.14[) 
and (|4.15[) . As a special case, this solution contains the 
discrimination of N symmetric linearly independent pure 
states with equal mutual overlaps, see Eqs. (14. 18|) and 
(I4.19j) . Moreover, we also obtained the solution for the 
discrimination of N symmetric equiprobable mixed qubit 
states, given by Eqs. (|4.21[) and (I4.22|) . and of N spe- 
cial symmetric mixed states of rank D spanning a joint 
Hilbert space of dimension ND, see Eq. (|4.28[) . 

(iii) We solved the optimization problem for a case of 
mixed qudit states that are complete in the sense that a 
weighted sum of their density operators is equal to the 
identity operator, and we also specified the general opti- 
mality conditions for a certain kind of partially symmet- 
ric states. The treatment of the optimization problem 
resulting for the latter case is left for further investiga- 
tions. 

In the final part of the paper we showed that there ex- 
ists a general relation between the solutions for optimum 
state discrimination in the two different discrimination 
strategies where either the rate Q of inconclusive results, 
or the overall error rate P e , has a fixed value. This rela- 
tion, expressed by Eqs. (I6.ip - (|6.3[) . holds for an arbi- 
trary number N of mixed states. It implies that by solv- 
ing the optimization problem in one of the two strategies, 
one can also obtain the solution in the other strategy. As 
an illustration we presented an example where for two 
mixed qubit states the maximum rate of correct results 
with a fixed error rate P e is derived from the solution 
for optimum state discrimination with a fixed rate Q of 
inconclusive results. 

In order to solve the optimization problem for a fixed 
probability Q of inconclusive results, we applied the op- 
erator conditions [25] determining the optimum measure- 
ment. As discussed in our paper, these optimality con- 
ditions, Eqs. (|2.6p and (|2.7[) . provide a very general 
approach for treating various optimized state discrimi- 
nation measurements, as far as only overall probabili- 
ties, averaged over all outcomes, are considered. In the 



appropriate limiting cases, they describe minimum-error 
discrimination, on the one hand, while on the other hand 
they refer to optimized maximum-confidence discrimina- 
tion provided that the maximum confidence is the same 
for each conclusive outcome, or to optimum unambigu- 
ous discrimination, respectively, in the special case when 
the maximum confidence is equal to unity. 

We note that related work has been done indepen- 
dently by E. Bagan and R. Munoz-Tapia (Barcelona) and 
G. A. Olivares-Renten'a and J. A. Bcrgou (New York) 

M. 
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Appendix A 

Here we present the derivation of the optimum mea- 
surement that discriminates two mixed qubit states when 
d = C 2 = C. Eq. d3U), valid for r - al > 0, yields 



T22 — Tn — 



a 2 - C(l - C) 
2a- 1 



, |ri 2 |=r n -a, (Ai) 



where we have to require that Tn > a and Tn > C > 
1/2 to ensure the positivity of the expressions under the 
square-root signs in Eqs. (|3.10[) - p,12j) . From Eqs. 
(|3.14j) and (|3.16|) we obtain the relations 

e iS = -e* £ 1+/ ? 2 = J^ = %^. (A2) 

21 11 — 1 |i 12 I 

Using Eq. (|A1|) in the last equality of Eq. (|A2I) we get 



1 2C-1 / 1 — 2|pi 2 



2 2 V 1 + 2|pi 2 | -2Q 



(A3) 



After calculating the matrix elements Tn and T 2 2 that 
result for this value of a, we obtain from Eqs. (|3 . 14[) and 
(13. 15)) the solutions 



y/(l - 2|pi 2 |)(l + 2\p 12 \ -2Q)± (gu - p 22 ) 
2(2(7-1) 



01/3 = 

(A4) 

Due to Eqs. (IA2|) and (|A4[) a solution with positive con- 
stants P\ and /?2 can only result when the two conditions 

Q<2\p 12 \, Q < P"P*-\pf ^ Qcr (A5) 
2 - IP12I 

are simultaneously fulfilled. In the second condition we 
used the relation (pn — P22) 2 = 1 — ^PuP22 and intro- 
duced the critical value Q cr of the failure probability. For 
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later purposes we note that 

2|pi2| < Qcr if |/9i2 1 < min{pn,p 2 2}, (A6) 
2|pi2 1 > Qcr if I/O12 1 > min{pn, / 922}- (A7) 

To see this, let us assume that \pi2\ < Pn < \, which 
leads to (~ - |/?i2 1) 2 > (| - Pn) 2 = | -P11P22. Hence it 
follows that P11P22 - IP12I 2 > I/O12 1 - 21/012 1 2 and conse- 
quently 2|pi 2 | < Q C r- Analogous considerations hold for 
the other possible cases. 

When the two conditions in Eq. (|A5[) are met, we 
obtain from Eq. (|3.8p the optimum detection operators 



n^Ap- 1 / 2 !^)^^- 1 / 2 (* = 1,2), (A8) 

where the states result from Eqs. (|3.11|> and (|3.12l) . 



> = (Vfc+Ih) + e-^v^ll^)) (A9) 



^2 



2C- 1 



with 6 = 



(V&^Th) + + (A10) 

1-Q 



V(l-2|p 12 |)(l + 2|p 12 |-2Q) 



(AH) 



For completeness, we also give an explicit expression for 
the failure operator. Because of Eqs. Q3.7P and (|3.10[) 
the latter can be written as 



n = Qp- 1/2 \Tt )(K \p 



-1/2 



Q 



\ir }(n \, (A12) 



where |7To) = \JQcrP 1 ^ 2 | 7r o) w ith (7ro|7ro) = 1. Here we 
took into account that 



TO 



PI 



V2 



and (tto|p Vo> = 75—, 

W cr 



(A13) 

where in the second equality use has been made of 
the matrix elements of p^ 1 with respect to the basis 

{h» 2 >}- 

If |/0i2 1 < mm{pu, P22} and therefore 2|pi 2 | < Q cr , 
the condition Q < 2 1 P12 1 is sufficient to guarantee that 
both inequalities in Eq. (|A5[) hold, which means that the 
detection operators given by Eqs. ((All) - (|A12|t describe 
the optimum measurement. Using Eq. (|2.12p together 
with the expressions for p\ and p 2 that ensue from Eqs. 
(|3.ip and (|3.2p when C\ = C2 = C, we obtain with the 
help of Eq. (|A8P the maximum probability of correct 

results pj ' given in Eq. (|3.23p . After combining this 
result with Eq. (f3~2TT|) we arrive at Eq. (HO!?)) . 

We now focus on the case that \p\ 2 \ > minjpn, P22}, 
where due to Eqs. (|A5|) and (j A7|) the detection opera- 
tors given by Eqs. (| A8|) - (|A12[) describe the optimum 
measurement only as long as Q < Q cr - In order to de- 
rive the solution for Q > Q cr , we assume without lack of 
generality that |p 12 | > pn, which implies that pn < p 22 
and iliai /',""" = C(l - Q) if Q > Q x , see Eqs. (IXT71) 
and (|3.20p . Equation (IA4|) reveals that in this case the 



parameter f3\ vanishes for Q = Q cr , which means that 
n! = 0. Therefore if Q cr < Q < Qi we have to search 
for a projective measurement with 



n 1 = o, n 2 = |7r 2 )(7r 2 |, n = /-n 2 



(A14) 



that maximizes P c at the fixed value Q and yields 
P™ ax (Qi) — C(l — Qi). Instead of using again the opti- 
mality conditions, we can solve the optimization problem 
directly, due to the simple structure of the detection op- 
erators. For this purpose we make the convenient general 
ansatz 



^2 



7ip- 1/2 h> 



2 lx V72>~ 1/2 h) 



(A15) 



where the phase x an d the nonnegative parameters 71 
and 72 have to be determined. The normalization condi- 
tion (7r 2 |7r 2 ) = 1 leads to the constraint 

71^22+72^11-2^7x72 \p 12 \ cos(x+0) = pnP22-|pi2| 2 = A, 

(A16) 

where we used the matrix elements of p 1 with respect 
to the basis {|^i), ^2)} and took into account that P12 = 
I/O12 le 1 ^. The fixed probability Q of inconclusive results 
can be expressed as 



l-Tr(pn 2 ) = l-(7i+ 72 ). 



(A17) 



From P c = 77 2 Tr(p 2 n 2 ) = Tr (p 2 p 1/2 U 2 p 1/2 ) we obtain 
the probability of correct results, 



P c = (l-C) 7l +C 72 , 



(A18) 



where Eq. (|3.2p has been used. In order to maximize 
P c we need to allow for the largest possible values of 71 
and 72 consistent with Eq. (| A16|) , that is we have to put 
cos(x + 4>) = 1, which means that \ = — 4>- Using Eqs. 
(|AT6| and (fAl"7) we then find that 



72 



(IP12IV71 + v^^/pi 



7i 



= l-Q-71. (A19) 



The equation arising from the second equality sign in 
Eq. (| A19|) can be solved to yield an explicit expression 
for 71 (Q). Making use of Eq. (|A18|) together with Eq. 
(|A17|) we get the solution 

prnax = p(l) = _ Q} _ ^ _ ^(Q), ( A20 ) 



which is explicitly given in Eq. (|3.26l) . Clearly, our 
derivation requires that 71 (Q) > 0. According to Eq. 
(|A19p the boundary case 71 = implies that A/pn = 
1 — Q and therefore Q = Qi, see Eq. (|3.17p . Hence Eq. 

(|A20)) is valid for Q cr < Q < Qi, yi eldin g P c (1) (Qi) = 
C(l — Qi), as required. Taking Eq. (|3.20p into account 
we arrive at the final result, Eq. (|3.24p . A direct calcu- 
lation shows that indeed Pc°\Qcr) = Pc^iQcr)- 

We still consider the question under which condition 
the relation C\ — C 2 = C holds true when the mixed 
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qubit states are represented with the help of the gen- 
eral ansatz p 1/2 = Pl/ 2 |l/'l/2>(^l/2l + 1 ^ P 2 2 !> See Ec l- 

(|3.27[) . After calculating the spectral representations of 
Pj = p^^-t^rjjPjp^ 1 / 2 (j = 1,2) and comparing the re- 



sult with Eqs. (|3.1I) and (13.21) we arrive at the condition 
(t?iPi) 2 — (j)2P2) 2 =1)1 — 1)2, leading to the explicit ex- 
pression 



(2C-1) 2 = 1 -- 



2V(l-^ 2 )(l-fI) 



1 + ^{l-pl)(l-pl)+p lP2 (l - 2\S\ 2 ) 

(A21) 

with 5 1 = (ipilipa). Taking into account that Ti{pip 2 ) = 
Tr {ppipfc)/ (V1V2) we find with the help of Eqs. (|3.27l) . 
(l3~T1) and (10) that 



l- PlP2 (l-2\S\ 2 ) _ C(l-C) (p? 1 +p| 2 -2| Pl2 | 2 ) + | Pl2 | : 



2 

(A22) 

Together with Eq. ()3.6|) and with the relation 

(77^ 1 — l) 2 = following from the condition for 

Ci = C 2 = C, Eqs. (|MT1) and (|Al2l) form a system 
of four equations for determining the values of i~)i , pi , P2 
and j S\ when C and the matrix elements pij are given. 



which arises from the positivity constraint IIo > with 
|^2| > \bi\- The resulting failure probability is given by 

Q=Tr(pno) = (l-iV|6 1 | 2 )mn- (B5) 

Taking into account that [Z, p] = it becomes obvious 
from Eqs. (|B2[) and (|B4|) that the equality in the first 
optimality condition, Eq. (|4.5|) . holds true if 



Z = «i Y]h)(n| + Z2 V \n)(n\ anda = -. (B6) 

1=1 !=m+l 

Now we turn to the equality in the second optimality 
condition, Eq. P~o]) . leading to — = 0, 

which implies that N(ri\Z\v) = (ri\pi\v) for I = 1, . . . ,d. 
After expressing pi with the help of Eqs. (|B1|) and (14. 13|) 
we arrive at the condition 

z i = ^ [™ c i fo i + ( d " m)c|6 2 ] + (i = 1, 2). 

(B7) 

Since Z is Hermitian we choose the phase of hi to coincide 
with the phase of Cj (i = 1,2). Equations (|B4[) and (|B5[) 
then lead to 



Appendix B 

We consider the discrimination of N symmetric qudit 
states pj (j = 1, . . . , N) described by 

pj=p\^m\ + ^ l i with \i )j ) = v«- i ^ 1 ), 

(Bl) 

where £ p < 1 and = ciJ^Zi h) + c 2 Ef=m+i h)> 
see Eq. (|4.13j) . Here / is the identity operator in the d- 
dimensional Hilbcrt space 7^ spanned by the N pure 
states \ipj). Using the properties of the symmetry op- 
erator V, given after Eq. ()4.2j) . we can determine the 



spectral representation of p = -k Pj, arriving at 



P = ri22\n)(ri\ +r 2 ^ \n){ri\, n=p\ 



l=l 



1-p 



(B2) 



with i — 1,2. The structure of pi suggests the ansatz 

m d 

Hi = l€f><«| with \v) = hJ2\n)+b2 X)|n>- (B3) 



1=1 



i=m+l 



To be specific, we assume that < |&i| < \b 2 \. We shall 
use the optimality conditions, Eqs. (|4.5|) and (|4.6|) . in or- 
der to investigate whether this ansatz yields the optimum 
solution for the case that Z — ap > 0. 

The equality in Eq. (|4.5|) implies that IIo cannot span 
the whole d-dimensional Hilbert space if Z — ap > 0. 
Using Eqs. (|4.7I) and (|B3I) this leads to the requirement 



N\b 2 \ 2 



1, U = (l-N\b 1 \ 2 )J2\r l )(r l \, (B4) 



1 c 2 , w ci /mri - Q 

2 = — 7=] — r> 0l = — ==1 — r with W = 



N \c 2 \' 



N ci 



From Eq. (|B7|) we obtain the relation 



l-p p\cic 2 \ ( |ci| d-m 

Z\ 7T^ = 77 I 171 



(B8) 



dTV N \ \c 2 \ w 



|ci| ( 1-P 
22 



diV 
(B9) 

yielding explicit expressions for 21 and z 2 . Using Eq. 
(|4.8[) we get the maximum probability of correct results, 
pma*\ Q = (v\pi\v) = P c (0) , where because of Eq. (|B8j 

P i 0) = ^ M C l| W + ( d " m )|c 2 |] 2 + ("^ + 

(BIO) 

In order to decide whether this is indeed the optimum so- 
lution we have to check whether the positivity constraints 
in Eqs. (|4"3j) and (|4"^|) are fulfilled. 

We start with the positivity constraint in Eq. 



which takes the form Z — — -fj§I > 0. Clearly, 

the constraint is satisfied provided that Nd(ipi\Z\ipi) > 
1 — p + pd. After inserting the explicit representation of 
the operator Z, this inequality can be transformed into 



m 2 |ci \ 4 +m(d-^m) \c\c 2 \ 



w\c 2 \ 



w\c 2 



-{d—m 



l 2 |c 2 | 4 >l. 



Due to the normalization condition in Eq. (|4.13p and to 
the fact that x + x^ 1 > 2 for any x it is obvious that Eq. 
(jBllj) always holds true, which implies that the positivity 
constraint in Eq. (|4.6[) is satisfied. 

The positivity constraint given in Eq. (I4.5[) . on the 
other hand, which results from the fact that Q is fixed 
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and different from zero, imposes the relevant restriction. 
Because of Eq. (|B6I) this constraint takes the form 



-2 



ar 2 = z 2 - z 1 — > 0, 

n 



(B12) 



where Zi and z 2 follow from Eq. (|B9I) . It turns out that 
the constraint is satisfied in the following general cases: 
(i) When p = 1 and therefore r, = |c,| 2 (i = 1,2), Eq. 
(fBT2"l) is fulfilled if w > |c 2 |/|ci|, which with the help of 
Eqs. (|B8|) and (14.131) leads to the condition 



<m(| Cl | 2 -|c 2 | 2 Hl-d|c 2 | 2 = 



(B13) 



(ii) When d = 2, and consequently m = d — m = 1, a 
straightforward calculation shows that Eq. (|B12[) holds 
true provided that w 2 > r 2 /r\ . Because of Eq. (|B8|) this 
yields the requirement Q < n — r 2 , or, using Eq. (|B2I) . 



<p(| Cl | 2 -|c 2 | 2 )=p(l~2| C2 | 2 ) 



(B14) 



(iii) When |ci| 2 = |c 2 | 2 = 4 for arbitrary values of p and 
d, implying that also n = r 2 = 4, Eq. (|B12I) is fulfilled 
if w = 1, that is, if Q = 0. This example belongs to 
the special case where the density operators resolve the 
identity operator, as discussed in Sec. V A. 



The preceding calculations have been performed as- 
suming that \b 2 \ > \bi\. Since Q has to be positive, 
we conclude from Eqs. (IB13[) and (|B14I) that this as- 
sumption is justified if \c 2 \ < |ci|. In the opposite case, 
l c 2| > |ci|, the derivation could be performed in a com- 
pletely analogous way, where the results arise from the 
previous ones when the indexes 1 and 2 are interchanged 
and when to is replaced by d — to, and vice versa. 
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